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Absfwct 

Helicopter-towed  electromagnetic  induction  sounding  systems  hove  recently 
been  used  for  the  remote  measurement  of  sea-ice  thickness  and  shallow  sea 
bathymetry.  An  inherent  aspect  of  this  sensing  technology  is  the  large  area,  the 
footprint,  over  which  an  individual  sounding  is  made.  This  report  gives  an 
overview  of  previous  footprint  size  assessments  and  presents  new  footprint  size 
estimates  determined  from  sea  ice  survey  data  and  analytical  calculations.  All 
footprint  determinations  were  found  to  be  in  reasonable  agreement.  For  a  vertical 
coaxial  coil  antenna  arrangement  the  apparent  footprint  diameter  was  found  to 
be  about  1 .25  times  the  antenna  height  above  the  conductive  surfoce,  and  for 
a  horizontal  coplanarcoil  configuration  the  ratio  is  about3.75times  the  antenna 
height. 


For  conversion  of  SI  metric  units  to  U.S./British  customory  units  of  measurement 
consult  ASTM  Standard  E380-89a,  Standard  PracticeforUso  of  the  International 
System  of  Units,  published  by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19103. 
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Footprint  Size  of  a 

Helicopten-Bome  Electromagnetic  Induction 
Sounding  System  Versus  Antenna  Altitude 

AUSTIN  KOVACS,  J.  SCOTT  HOLLADAY  AND  CLYDE  J.  BERGERON,  JR. 


INTRODUCTION 

Ice  covers  large  areas  of  the  Arctic  and  Antarctic 
seas.  The  thickness  and  distribution  of  this  veneer 
are  of  considerable  importance  for  efficient  ship 
routing,  design  of  offshore  structures,  and  scientific 
sea  ice  dynamic-thermodynamic  studies.  Before 
1985,  many  techniques  were  evaluated  for  the  re¬ 
mote  measvirement  of  sea-ice  thickness.  Most  met 
with  limited  success  because  of  losses  associated  with 
the  higj\  brine  inclusicm  content  and  related  high  con¬ 
ductivity  of  sea  ice. 

In  1985,  a  conventional  helicopter-towed  geophysi¬ 
cal  four-frequency  electromagnetic  (EM)  induction 
soimding  system  was  first  used  to  estimate  sea-ice 
thickness  (Kovacs  et  al.  1987).  The  EM  soimding  sys¬ 
tem  works  according  to  the  principles  of  electromag¬ 
netic  induction,  where  the  dectromagnetic  fields  are 
governed  by  diffusion  radier  than  a  wave  equation. 
A  transmitter  coil  is  excited  by  a  sinusoidal  electric 
current.  This  produces  an  alternating  magnetic  field 
in  space  so  that  electromagnetic  currents  are  induced 
in  any  nearby  conductors  (e.g.,  seawater).  These  in¬ 
duced  currents  in  turn  produce  a  secondary  mag¬ 
netic  field  that  is  sensed  by  the  receiver  coil. 

Located  halfway  between  the  transmit  and  the  re¬ 
ceiver  antenna  coils  is  a  "bucking  coil"  At  this  loca¬ 
tion  the  primary  field  is  about  eight  times  stronger 
than  at  the  receiver  coil,  but  the  secondary  field  is 
essentially  the  same.  Therefore,  the  bucking  coil  is 
designed  to  have  eight  times  less  effective  area  than 
the  receiver  coil  so  that  simple  subtraction  of  the  re¬ 
ceived  signal  at  the  bucking  coil  from  the  received 
signal  at  the  receiver  coil  "bucks  out"  the  primary 
field's,  leaving  a  slightly  attenuated  secondary  com¬ 
ponent.  Using  this  technique,  tihe  secondary  field  can 
be  measured  with  high  precision;  its  in-phase  and 
quadrature  components  are  recorded  digitally  in 


parts  per  million  (ppm)  of  the  primary  magnetic  field 
with  a  precision  of  about  1  ppm.  These  data  are  used 
to  determine  the  hei^t  of  the  antenna  housing  (bird) 
above  the  seawater  surface  and  its  conductivity.  At 
the  same  time  the  bird's  altitude  above  the  ice  sur¬ 
face  is  measured  by  a  laser  altimeter.  Thus  the  ap¬ 
parent  pack-ice  tiiickness  is  obtained  simply  from 
the  difference  in  these  two  distances. 

The  1985  EM  sounding  field  trials,  made  over  the 
winter  Beaufort  Sea  ice,  were  encouraging.  Further 
EM  system  refinements  led  to  halving  the  towed  sen¬ 
sor  platform  length  from  about  6.5  to  3.5  m  (Kovacs 
and  HoUaday  1990)  and  then  to  the  development  in 
1990  of  a  wideband  system  having  a  firequency  range 
from  about  5  to  200  kHz  and  the  unique  capability 
to  process  the  data,  using  a  one-dimensional  layered- 
halfspace  inversion  routine,  and  display  the  ice  thick¬ 
ness  in  real  time  at  >10  Hz.  At  the  inception  of  tins 
EM  test  and  evaluation  program,  a  persistent  ques¬ 
tion  was:  How  large  was  the  area  over  which  an  in¬ 
tegrated  ice  thickness  measurement  was  made?  For 
the  purpose  of  this  discussion,  the  basic  unit  of  reso¬ 
lution  will  be  called  the  footprint.  This  wUl  be  loosely 
defined  as  die  horizontal  distance  along  a  survey  line 
from  which  90%  of  the  groimd  response  of  the  sys- 
ton  arises.  General  rule-of-thumb  estimates  that  were 
available  in  the  geophysical  community  in  1985  gave 
the  footprint  of  a  helicopter-towed  electromagnetic 
(HEM)  induction  sounding  system  as  being  about  2 
to  3  times  the  antenna  height  above,  for  example,  a 
conductive  seawater  surface.  The  size  of  the  foot¬ 
print  is  important  to  understanding  the  degree  to 
which  HEM  sounding  spatially  smoothes  out  the  un¬ 
dulating  seawater  relief  associated  with  ice  thickness 
variations.  In  short,  a  large  footprint  will  tend  to 
smooth  or  subdue  the  thickness  of  pressure  ridge 
keels  and  ice  relief  that  has  a  width  on  the  order  of 
the  sounding  system's  footprint  or  smaller.  In  this 


report  we  present  recent  assessments  of  the  HEM 
footprint  size. 


TRANSMITTED  HELD  SHAPE 

A  HEM  system's  antenna  housing  (bird)  is  gener¬ 
ally  suspended  about  30  to  40  m  below  the  helicop¬ 
ter  and  is  towed  about  20  to  30  m  above  the  surface. 
Inside  the  dgar-shaped  bird  are  circular  coil  pairs  (a 
transmit  Tx  coil  and  a  receiver  Rx  coil).  The  coil  pairs 
are  generally  arranged  as  shown  in  Figure  1,  with  a 
bucking  coil  located  in  the  same  orientation  halfway 
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Figure  1.  Three  typical  HEM  system  coil  arrange¬ 
ments.  Arrows  indicate  the  magnetic  dipole  orien¬ 
tation  (from  Kovacs  et  al.  1987). 
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Figure  2.  Primary  electromagnetic  field  direction  for  two  transmitters. 
Arroivs  indicate  the  magnetic  dipole  antenna  orientations  (from  Kovacs  et 
al.  1987). 


between  them.  The  direction  (but  not  the  magnitude) 
of  the  magnetic  field  lines  generated  by  the  coaxial 
and  horizontal  coplanar  transmitter  is  ^own  in  Fig¬ 
ure  2.  In  these  diagrams  the  bird  axis  is  horizcattal. 
The  vertical  coplanar  coil  configuration  generates  a 
field  distribution  identical  to  that  of  the  coaxial  coil 
but  rotated  90°  into  the  page.  A  more  important  as¬ 
pect  of  the  effect  of  coil  configuration  used  in  sea  ice 
sounding  is  the  pattern  of  the  currents  induced  in 
the  seawater  by  the  Tx  coil,  because  it  is  these  in¬ 
duced  currents  fiiat  generate  the  signal  measured 
and  analyzed  by  the  HEM  system.  The  zone  in  which 
the  strongest  current  flow  occurs  will  be  the  region 
that  contributes  most  to  the  response  measured  by 
the  HEM  system.  Therefore,  the  smaller  the  area  of 
high  current  flow  the  smaller  the  footprint. 

The  model  used  by  Liu  (1989)  to  calculate  the  sur¬ 
face  current  distributions  shown  in  Figure  3  was  a 
highly  conducting  sheet  whose  effective  depth  de¬ 
pends  on  the  frequency.  The  angles  of  the  arrows 
indicate  the  directions  of  the  currents,  while  their 
length  denotes  fire  relative  strength  of  the  current  at 
that  position.  At  high  frequencies  (>10  kHz),  these 
currents  lie  close  to  the  seawater  surface  and,  since 
their  distance  from  the  antenna  coil  can  be  estimated 
accurately,  the  use  of  high  frequencies 
permits  an  accurate  estimate  of  the  an¬ 
tenna  altitude  above  the  water  surface. 
At  lower  frequencies,  this  pattern  be¬ 
comes  more  diffuse  and  is  centered  at 
a  greater  depth  beneath  the  water  sur¬ 
face.  This  is  why  there  is  better  (deeper) 
bathymetric  sampling  by  a  system  op¬ 
erating  at  lower  frequencies,  e.g.,  at  50 
Hz.  A  strong  concentration  of  current 
located  immediately  below  the  trans¬ 
mitter  (at  X  and  y  =  0)  is  shown  in  Fig¬ 
ure  3a.  Here  the  zone  of  maximum  cur¬ 
rent  strength  lies  in  an  ellipse  with  its 
short  axis  lying  along  the  bird  axis  (x- 
direction)  and  its  long  axis  perpendicu¬ 
lar  to  the  bird  axis  (y-direction).  Two 
zones  of  weak  response  precede  and 
follow  the  bird  along  its  axis.  The  bird 
height  assumed  for  this  diagram  and 
the  next  was  30  m. 

The  induced  surface  currents,  as 
calculated  by  Liu  (1989),  for  the  hori¬ 
zontal  coplanar  transmitter  mode  are 
shown  in  Figure  3b.  In  this  configura¬ 
tion,  the  currents  are  strongest  in  a  ring 
with  a  radius  of  about  two-thirds  the 
flight  height.  Directly  below  the  trans¬ 
mitter  there  is  a  narrow  zone  of  zero 
response.  The  response  dies  off  sharply 
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Figure  3.  Idealized  surfacial  electric  current  pattern  (from  Uu  1989). 
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Figure  4.  Representative  surfacial  current  density  profile  (from  Uu 
1989). 


as  the  distance  from  the  transmitter  becomes 
larger  than  about  one  flight  height. 

The  current  distribution  for  Ihe  vertical  co- 
planar  configuration  can  be  obtained  by  tiun- 
ing  Figure  3a  90°  clockwise  and  considering 
the  flight  direction  to  lie  along  the  new  r-axis 
(i.e.,  off  to  the  right).  The  elliptical  zone  of 
rrwximum  response  would  now  have  its  long 
axis  along  the  flight  direction  and  its  short 
axis  perpendicular  to  this.  This  mearxs  that 
the  vertical  coplanar  configuration  tends  to 
average  features  encountered  along  the  flight 
line  over  a  considerably  longer  distance  than 
does  the  coaxial  configuration,  and  that  it  is 
relatively  insensitive  to  features  that  lie  off  to 
the  side  of  the  flight  line  by  more  than  about 
one  flight  height. 

The  current  density  in  profile  from  along 
AA'  and  BB'  in  Figure  3a  is  shown  in  Figure 
4.  The  dominant  response  along  AA'  in  Fig¬ 
ure  4a  is  seen  to  arise  from  a  zone  about  1.3 
times  the  flight  height  in  width  along  the 
flight  line.  Side  lobes  are  also  present  in  which 
the  currents  actually  run  coimter  to  their  di¬ 
rection  in  the  main  zone.  Slicing  die  distribu¬ 
tion  perpendicular  to  the  flight  line  along  BB' 
yields  the  pattern  shown  in  Figure  4b.  Here 
the  response  is  seen  to  arise  primarily  from  a 
zone  about  four  flight  heights  wide. 

The  vertical  coplanar  current  density  pro¬ 
files  can  be  obtained  by  simply  replacing  Fig¬ 
ure  4b  with  Figure  4a.  This  again  demon¬ 
strates  the  reduced  resolution  along  the  flight 
path  of  about  four  times  less  than  for  the  co¬ 
axial  case. 
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Figure  5.  Representative  surfacud  current  density  profile  along  AA' 
and  BB"  in  Figure  3b  (fromLiu  1989). 


The  current  density  profile  for  the  horizontal  copla- 
nar  case  is  shown  in  Figure  5.  Since  the  current  dis¬ 
tribution  is  axially  symmetric,  profiles  AA'  and  BB' 
are  identical.  The  null  point  below  the  transmitter 
shows  up  clearly,  as  does  the  peak  of  the  current 
ring  at  plus  or  minus  two-thirds  the  flight  height. 

HEM  FOOTPRINT  SIZE 

During  the  course  of  developing  and  evaluating 
HEM  technology  for  the  measurement  of  sea-ice 
thickness,  a  number  of  studies  were  made  to  better 
define  the  HEM  footprint  size.  An  empirical  assess¬ 
ment  was  made  by  Kovacs  et  al.  (1987),  in  which  a 
HEM  system  was  flown  above  a  grid  area  on  a 
multiyear  ice  floe.  Drillhole-measured  ice  thickness 
varied  from  about  2  to  5  m  within  the  grid  area.  The 


drillhole-measured  ice  thickness  was  then 
averaged  within  a  successively  enlarged 
rectangular  area  under  the  HEM  antenna 
until  the  averaged  drillhole-measured  ice 
thickness  agreed  with  the  value  estimated 
from  the  measured  HEM  resptmse.  This  as¬ 
sessment  indicated  that  for  the  coaxial  coil 
arrangement  the  apparent  footprint  diam¬ 
eter  is  about  1.25  times  the  bird's  height. 
Since  the  shape  of  the  induced  current  field 
is  definitely  not  rectilinear  and  the  number 
of  drillhole  measurements  was  limited,  this 
footprint-altitude  ratio  estimate  may  not  be 
very  accurate. 

At  about  the  same  time  a  numerical 
analysis  of  HEM  footprint  size  was  carried 
out  by  Liu  and  Becker  (Becker  et  al.  1987, 
Liu  1989,  Liu  and  Becker  1990).  In  this  assessment 
the  fraction  of  the  total  HEM  response,  generated  by 
the  induced  current  flowing  at  ^e  surface,  within  a 
square  box  centered  beneath  the  transmitting  an¬ 
tenna  was  estimated.  The  length  of  the  box's  sides 
was  varied  to  estimate  the  fractional  response  for 
different  "footprint"  sizes  for  a  transmitter  situated 
30  m  above  a  seawater  surface.  The  footprint  esti¬ 
mates  for  the  coaxial  and  horizontal  coplanar  modes 
are  shown  in  Figure  6  (Liu  1989).  At  the  90%  fraction 
contribution  the  apparent  footprint-altitude  ratio  is 
shown  to  be  1.35  for  the  coaxial  and  3.73  for  the  hori¬ 
zontal  coplanar  antenna  system.  The  former  is  in 
good  agreement  with  the  above  assessment  of 
Kovacs  et  al.  (1987).  Nevertheless,  since  the  induced 
current  system  is  definitely  not  rectilinear,  this  foot¬ 
print  estimator  may  also  overestimate  the  footprint- 
altitude  ratio  by  averaging  zones  of  high  current  den- 


Figure  6.  Estimated  snrfacial  atrreid 
density  contribution  at  a  30-m  eleva¬ 
tion  Jbr  a  coaxial  and  a  coplanar  HEM 
antenna  coil  orientation  versus  an  ide¬ 
alized  square  footprint  width  to  an¬ 
tenna  height  ratio.  Note  that  Jbr  the  co¬ 
planar  coil  orientation,  about  15%  of  the 
measured  response  can  be  expected  from 
a  box  having  a  width  equal  to  the  bird 
height  and  about  90%  can  be  expected 
0  1  2  3  4  5  from  a  box  having  a  width  of  3.7  times 

Square  Width  /Antenna  Height  Ratio,  w/h  the  bird  height  (after  Uu  1989). 
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Figure  7.  Measured  response  of  a 
coaxial  HEM  system  flown  (a) 
from  over  thick  first-year  sea  ice 
out  over  a  lead  mth  thin  ice  and 
(b)  from  the  thin  lead  ice  back  over 
the  thicker  first-year  sea  ice. 


sity  (and  therefore  contributing  to  a  stronger  system 
response)  with  zones  of  much  weaker  current.  Liu 
(1989)  did  not  calculate  an  apparent  footprint  size 
for  the  vertical  coplanar  system,  but  by  exammation 
of  the  currait  density  profiles  and  comparison  with 
the  existing  estimates,  a  reasonable  value  would  be 
approximately  four  times  the  anterma  height. 

The  wideband  HEM  induction  sounding  system 
was  field  tested  near  Prudhoe  Bay,  Alaska,  in  April- 
May  1990.  This  system  could  be  operated  in  either  a 
coa^dal  or  a  horizontal  coplanar  mode  at  a  sampling 
rate  of  50  Hz.  An  exampk  of  each  mode  will  be  dis¬ 
cussed  as  related  to  the  footprint-altitude  ratio.  The 
frequencies  used  (30  to  100  kHz)  were  high,  so  ti\at 
the  majority  of  (he  induced  current  flow  occurred 
near  the  surface  of  the  seawater. 

A  site  was  selected  where  a  step  in  sea-ice  thick¬ 
ness  occurred  from  fairly  thick  (1.^  m)  first-year  ice 
to  very  rtrin  (0.1  m)  lead  ice.  The  HEM-measured 


response  over  this  transition  was  plotted,  and  the 
start  and  end  of  the  HEM  response  arising  from  the 
step  change  in  ice  thickness  was  estimated.  This  was 
then  compared  to  the  average  height  of  the  HEM 
^stem  and  the  aircraft  speed  (65  knots)  to  estimate 
the  flight  distance  vs.  fiducial  number.  The  raw  HEM 
data  were  analyzed.  No  anti-alias  Bitering  or  spheric 
noise  removal  was  performed,  although  the  data 
were  baselined  for  drift  removal  and  calibrated  to 
read  in  parts  per  million.  The  apparent  footprint  for 
each  coil  configuration  was  then  estimated. 

The  response  of  the  coaxial  system  at  the  first- 
year/lead  ice  transition  for  two  overflights  is  shown 
in  Figure  7.  The  fiducial  width  of  the  transition  zone 
has  been  marked  (dashed  line),  as  has  the  central 
position,  which  represents  die  location  of  the  first- 
year  ice  edge  or  step  in  ice  thickness.  For  the  data  in 
Figures  7a  and  7b  the  analysis  indicates  an  apparent 
footprint  of  13.9  and  13.2  m,  respectively,  or  0.8  and 
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Figure  8.  Measured  response  of  a  coaxial  HEM  system  as  it  was  flown  from  aver 
thick  first-year  sea  ice  out  over  thin  lead  ice  and  then  back  aver first-year  sea  ice. 


b.  1991. 

Figure  9.  Ice  thickness  distribution  along  the  199G  and  1991  survey  lines. 
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0.9  times  the  antenna  elevation  above  the  seawater. 
These  values  cire  appreciably  less  than  those  deter¬ 
mined  by  Kovacs  et  al.  (1987)  of  1.25  and  Liu  (1989) 
of  1.35  for  a  coaxial  coil  system.  However,  these  low 
footprint-alhtude  ratios  may  be  off  by  up  to  20% 
due  to  an  uncertainty  in  the  GPS-estimated  flight 
speed  at  the  instant  of  passing  over  the  step  area. 

The  respcttise  of  a  horizontal  coplanar  coil  system 
flown  over  a  narrow  lead  is  shown  in  Figure  8.  At 
an  estimated  flight  speed  of  65  knots  the  footprint 
works  out  to  be  approximately  66  m  at  an  average 
antenna  elevation  of  18  m  above  the  seawater.  This 
results  in  an  apparent  footprint-altitude  height  ratio 
of  3.7,  which  is  in  agreement  with  the  cissessment  of 
Liu  (1989)  of  3.7  at  the  90%  response  contribution 
level  as  shown  in  Figure  6. 

During  the  1990  and  1991  HEM  sea-ice  thickness 
measurement  test  and  evaluation  program,  control 
lines  were  laid  out  across  the  sea  ice.  The  1990  line 
was  1.3  km  long  and  extended  across  first-year  sea 
ice  from  1.3-m-thick  ice  in  a  lead  to  11-m-thick  ice  in 
a  pressure  ridge.  The  1991  line  was  1.0  km  long  and 
crossed  both  first-year  and  multiyear  sea  ice  that  var¬ 
ied  from  about  1.6  to  14  m  in  thickness.  Along  each 
line  the  snow  cover  plus  ice  thickness,  hereafter  called 
the  ice  thickness,  was  determined  by  drillhole  mea¬ 
surements  made  at  5-m  intervals.  The  drillhole-mea¬ 
sured  ice  thickness  along  the  1990  and  1991  survey 
lines  is  shown  in  Figures  9.  The  mean  thickness  and 
standard  deviation  of  the  ice  along  the  1990  survey 
line  were  2.65  and  1.48  m,  respectively,  and  for  the 
1991  survey  line  they  were  3.02  and  1.85  m,  respec¬ 
tively. 

Two  HEM  ice-thickness  sounding  runs  were  made 
over  the  1990  survey  line  using  the  horizontal  copla¬ 
nar  coil  configuration.  The  average  height  of  thie  bird 
above  the  seawater  during  pass  1  was  20.7  m  and 


for  pass  2  it  was  22.8  m.  The  mean  and  standard 
deviation  for  pass  1  were  2.66  m  and  1.16  m,  respec¬ 
tively,  and  for  pass  2  they  were  2.72  and  1 .22  m.  There 
is  renaarkable  agreemait  between  the  mean  drillhole- 
measured  and  HEM-estimated  ice  thickness  values. 
This  agreement  was  a  result  of  the  short  HEM 
flightlines,  during  which  time  nonlinear  system  drift 
was  nrinimal,  and  the  resulting  HEM  data  were 
therefore  of  the  highest  quality.  An  example  of  the 
HEM-estimated  ice  thickness  for  pass  1  vs.  the 
drillhole-measured  values  is  shown  in  Figure  10. 
Note  that  the  HEM  data  do  not  fully  reveal  the  depth 
of  the  pressure  ridge  keels.  This  is  a  result  of  the 
spatial  smoothing  of  ice  relief,  which  occurs  as  a  re¬ 
sult  of  the  wide  HEM  footprint. 

An  assessment  of  the  horizontal  coplanar  coil  foot¬ 
print  size  was  made  as  follows.  The  drillhole-mea¬ 
sured  and  the  HEM-estimated  ice  thickness  vs.  dis¬ 
tance  along  the  1990  survey  line  were  analyzed  to 
reveal  the  normalized  autocorrelation  function  of  ice 
thickness  vs.  the  lag  distance  as  outlined  in  Buznev 
and  Dubovtsev  (1971)  and  Rothrock  (1986).  Of  in¬ 
terest  is  the  lag  distance  at  which  the  autocorrelation 
function  curve  has  its  first  zero  crossing.  This  repre¬ 
sents  the  minimum  distance  between  statistically  in¬ 
dependent  ice  thickness  measurements.  For  the 
drillhole  measurements  and  the  pass  1  and  pass  2 
HEM  data  this  distance  is  about  60,  70,  and  80  m, 
respecti veJy,  as  shown  in  Figure  1 1 .  A  running  aver¬ 
age  of  the  drillhole  data  was  then  rr  ade  in  which  the 
data  bin  size  was  progressively  increased  until  the 
autocorrelation  function  for  the  running  average  data 
set  passed  through  zero  at  the  same  lag  distance  as 
the  HEM  data. 

The  final  bin  size  resulting  from  this  process  pro¬ 
vided  the  number  of  drillhole  measurements  and  the 
linear  distance  along  the  survey  line  over  which  the 


Figure  10.  HEM-determined  vs.  drill¬ 
hole-measured  ice  thickness  along  the 
1990  survey  line. 
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Figure  11.  Auto-correlation  function  cf  ice  thickness  vs.  lag  distance  for  (a)  the  drillhole- 
measured  and  (b  and  c)  the  HEM-determined  ice  thidmess  along  the  1990  survey  line. 


measurements  were  averaged  to  obtain  a  zero  lag 
distance  comparable  to  the  HEM  data.  For  passes  1 
and  2  d\e  distances  were  73  and  100  m,  resp^vely. 
Since  the  average  antenna  height  above  the  seawa¬ 
ter  for  passes  1  and  2  was  20.7  and  2Z8  m,  respec¬ 
tively,  the  apparent  footprint-altitude  ratio  for  each 
pass  was  3.5  and  4.0,  respectively. 

Three  HEM  ice-thiclmess  soimding  runs  were 
made  over  the  1991  survey  line.  The  average  bird 


height  during  passes  1, 2,  and  3  was  23.0, 22.1,  and 
24.0  m,  respectively.  The  mean  and  standard  devia¬ 
tion  for  the  HEM-estimated  ice  thickness  were:  for 
pass  1, 2.90  and  0.54  m;  for  pass  2,  Z90  and  0.59  m; 
and  for  pass  3, 3.01  and  1.85  m.  The  agreement  be¬ 
tween  these  mean  ice-thickness  values  and  those 
determined  from  the  drillhole-measured  data  of  3.02 
m  is  very  good,  for  the  same  reascoi  previously  given 
for  the  good  comparison  between  the  1990  data  sets. 
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Figure  11  (cant’d). 


An  example  of  the  HEM-determined  ice  thickness 
for  pass  2  vs.  the  drillhole-measured  values  is  shown 
in  Figure  IZ  While  the  mean  thickness  as  determined 
by  the  HEM  system  (2.90  m)  and  the  drillhole-mea¬ 
sured  value  (3.02  m)  are  in  good  agreement,  it  is  ap¬ 
parent  from  a  comparison  of  Figures  10  and  12  that 
more  smoothing  occurred  in  the  1991  HEM  profile 
data.  This  is  partly  because  the  a\  erage  height  of  the 
bird  above  Ae  ice/ water  interface  was  2.6  m  more 
in  1991  than  in  1990.  The  additional  flight  height 
translates  into  a  slightly  larger  footprint  (Fig.  6)  and 
thus  more  spatial  smoothing  of  the  ice  relief.  The  ice 


thickness  variations  along  the  1991  survey  line  were 
also  narrower  and  occurred  at  a  high  hequency.  Over 
ice  of  uniform  thickness,  the  IffiM-estimated  ice 
thickness  is  not  dependent  on  bird  altitude  and  thus 
footprint  size.  However,  over  deformed  ice,  the 
HEM-estimated  thickness  of  an  individual  sea-ice 
pressure  ridge  is  quite  dependent  on  the  bird  alti¬ 
tude.  At  a  first-year  pressure  ridge  a  distortion  of 
the  EM  fields  occurs  due  to  fire  block  structure  and 
irregular  shape  of  the  ice  formation.  In  addition, 
where  the  ridge  is  narrower  than  the  footprint,  a  1- 
D  inversion  of  the  EM  data  does  not  provide  good 


Figure  12.  HEM-determined  vs.  drillhok- 
measured  ice  thickness  along  the  1991  sur¬ 
vey  line. 
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Figure  13.  Maximum  HEM-determined  pressure  ridge  ice  thick¬ 
ness  vs.  bird  elevation  above  the  surface. 


draft  estimates.  These  factors  affected  the 
HEM-estimated  ice  thickness  profiles,  but 
more  so  along  the  1991  survey  line  with  its 
narrower  ridges. 

The  effect  of  bird  elevadtm  and  therefore 
footprint  size  is  illustrated  in  Figure  13,  which 
shows  the  maximum  HEM-determined  ice 
thickness  over  tlw  thickest  ice  feature  shown 
in  Figure  9a  vs.  the  bird  height  above  the  sur¬ 
face  for  eight  soimding  passes  with  the  co- 
planar  antenna  orientation.  Except  for  the 
outlier,  at  the  22.3-m  altitude,  the  data  show 
the  decreasing  maximum  ice  thickness  trend 
vs.  antenna  height  relationship  that  would 
be  expected  from  the  footprint  effect.  The 
analysis  for  the  regression  line  shown  pass¬ 
ing  through  the  data  did  not  include  the  out¬ 
lier,  as  this  data  point  was  collected  slightly 
off  to  one  side  of  the  survey  line.  The  slope  of 
the  line  in  Figure  13  is  0.40  ±  0.12.  This  value 
is  in  agreement  with  the  value  of  0.3  ttiat  can 
be  estimated  from  Liu's  (1989)  footprint  analysis. 

As  with  the  1990  data,  a  normalized  autocorrela¬ 
tion  analysis  was  run  on  the  1991  data.  The  results 
are  shown  in  Figure  14  for  the  drillhole  data  and 
passes  1, 2,  and  3.  Here  the  minimum  distances  be¬ 
tween  statistically  independent  samples,  as  de¬ 
fined  by  where  the  autocorrelation  function  passes 
through  zero,  are  about  35, 132, 129,  and  119  m  for 
the  drillhole  measurements  and  passes  1, 2,  and  3, 
respectively.  The  running  average  autocorrelation 
analysis  of  the  drillhole  data,  as  previously  described, 
indicated  that  zero  lag  distances  or  "footprint" 
lengths  of  90,  88,  and  70  m,  respectively,  matched 
the  zero  lag  distance  for  passes  1, 2,  and  3.  The  aver¬ 
age  horizontcil  coplancir  coil  antenna  heights  above 
the  seawater  for  passes  1,  2,  and  3  were  23, 22,  and 
23.5  m,  respectively.  Therefore,  the  apparent  foot- 
print-altitude  ratio  for  each  pass  was  3.9,  4  0,  and 
3.4,  respectively. 

The  variations  in  the  estimated  footprint-altitude 
ratios  may  be  partly  due  to  bird  flight  path  devia¬ 
tions  away  from  the  survey  line  and  to  the  use  of  a 
single  line  of  drillhole  ice-thickness  measurements 
to  analyze  HEM  data  obtained  from  the  integration 
of  the  ice  thickness  over  a  much  wider  swath.  Nev¬ 
ertheless,  the  average  of  the  1990  and  1991  footprint- 
altitude  ratios  is  3.8.  This  value  is  in  very  good  agree¬ 
ment  with  the  3.7  ratio  determined  by  Liu  (1989)  and 
can  be  further  verified  as  follows. 

Bergeron  and  Michel  (1986)  calculate  the  footprint 
of  a  horizontal  coplanar  coil  pair  HEM  induction  sys¬ 
tem  in  the  context  of  the  Modified  Image  Method 
(MIM)  representation  of  the  secondary  field.  At  50% 


of  the  total  signal,  the  footprint  diameter  was  nu¬ 
merically  determined  to  be  twice  the  altitude.  Re¬ 
cently  the  functional  form  of  a  2-D  HEM  impulse 
function  for  the  horizontal  coplanar  coil  configura¬ 
tion  was  published  (Bergeron  et  al.  1992).  The  model 
for  the  impulse  function  is  fl\e  secondary  field  pro¬ 
duced  by  a  line  source  current  induced  in  the  image 
plane  by  the  primary  HEM  field.  This  is  the  same 
model  used  by  Liu  (1989)  except  tiiat  the  surface  cur¬ 
rent  sheet  is  placed  one  complex  screening  length 
below  the  conducting  surface. 

The  expressions  for  the  normalized  line  impulse 
function  of  the  MIM  representation  of  the  second¬ 
ary  field  for  the  horizontal  coplarwr  (HCP)  and  ver¬ 
tical  coaxial  (VCX)  coil  configurations  are  as  follows: 

HCP. 

(l  +  rf  ^ 

VCX: 

y  _  l-1.5r^ 

(i-f" 

where 

r  =-^ 
h* 

and 

h*  =  liait  +  8'ff 

in  which  x  is  the  horizontal  coordinate  of  the  line 
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source,  oriented  in  the  y  direction,  is  the  altitude 
of  the  coincident  source-receiver  coUs  above  the  con¬ 
duction  structure,  and  Seff  is  the  complex  effective 
skin  depth  of  the  conducting  structure,  which  makes 
h*  the  complex  altitude.  Because  diese  expressions 
contain  both  the  altitude  and  die  effective  sl^  depth, 
the  impulse  response  can  be  range-dependent. 

From  their  results  we  now  calculate  die  width  of 
t*"  line  impulse  function  at  10%  of  the  mtodmum 


value  or  90%  of  total  signal  level  for  die  HCP  coil 
configuration  and  find  it  to  be  3.7  h  where  h  is  the 
bird  altitude.  This  width  is  a  measure  of  the  2-D 
HEM  footprint.  For  the  VCX  coil  arrangement,  die 
calculation  of  the  10%  line  width  of  the  90%  2-D 
footprint  of  the  HEM  impulse  function  gives  a  value 
of  1.32  h.  These  results  agree  widi  Liu's  (1984)  and 
with  those  derived  from  field  data  as  presented  in 
this  report. 


b.  Pass  1. 

Figure  14.  Auto-correlation  function  of  ice  thickness  vs.  lag  distance  for  (a)  the  drillhole- 
measured  and  (b  through  d)  the  HEM^etermined  ice  thickness  along  the  1991  survey  line. 
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Figure  14  (cant'd).  Auto-correlation  function  cf  ice  thickness  vs.  lag  distance  for  (a)  the 
drillhole-measured  and  (b  through  d)  the  HEM-determined  ice  thickness  along  the  1991 
survey  line. 


CONCLUSIONS 

In  this  report  assessments  are  given  for  the  appar¬ 
ent  footprint-altitude  ratio  for  coaxial  and  horizon¬ 
tal  coplanar  coil  HEM  sounding  systems.  The  find¬ 
ings  indicate  that  the  apparent  footprint  diameter  of 
a  coaxial  and  a  horizontal  coplanar  HEM  system  is 
respectively  about  1.25  and  3.75  times  the  antenna 


height  above  the  sea  surface  and  that  the  smallest 
HEM  footprint  is  obtained  with  the  use  of  a  coaxial 
coU  arrangement.  In  principal,  this  configuration 
would  allow  for  a  better  assessment  of  ridge  keel 
depth  and  geometry,  but  not  ice  volume.  However, 
the  reason  for  not  using  a  coaxial  coil  configuration 
is  that  the  signal-to-noise  ratio  becomes  less  favor¬ 
able  by  a  factor  of  about  four. 
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